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^ ; Abstract 



Two-parametric eddy viscosity (TPEV) and other spectral characteristics 
of two-dimensional (2D) turbulence in the energy transfer sub-range are cal- 
culated from direct numerical simulation (DNS) with 512 2 resolution. The 
DNS-based TPEV is compared with those calculated from the test field model 
(TFM) and from the renormalization group (RG) theory. Very good agree- 
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ment between all three results is observed. 
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Two-dimensional incompressible turbulent flows are described by the vorticity equation: 

oc 9(no_ 2 

m + d(x, y ) ~ PoVC (1) 

where ( is fluid vorticity and vq is molecular viscosity. It is well known that the existence 
of inviscid invariants / d 2 x( 2n of (|1]) results in the flux of energy towards the largest spatial 
scales. The presence of this inverse cascade complicates the large-scale description of 2D 
flows and requires refinement of the classical hydrodynamic notion of "eddy viscosity." The 
concept of eddy viscosity is well defined for 3D turbulent flows, where energy cascades 
towards the smallest flow scales where it is dissipated. To achieve an adequate coarse- 
grained description of 3D flow, one can introduce increased "effective" dissipation at large 
scales which accounts for the unresolved dissipation. 

In 2D flows, the inverse flux of energy at large scales and enstrophy dissipation at small 
scales make the eddy viscosity concept more subtle. It was argued by Kraichnan |J that, 
in Fourier space, a 2D eddy viscosity should include two parameters: a cutoff wave number 
k c (which essentially determines the size of the coarse grain), and the wave number of a 
given mode, k. The two-parameter eddy viscosity (TPEV), denoted by u(k\k c ), describes 
the energy exchange between a given resolved vorticity mode with the wave number k and 
all subgrid, or unresolved, modes with k > k c ; it provides correct account for the energy and 
enstrophy fluxes between resolved and unresolved scales. The TPEV is derived from the 
evolution equation for the spectral enstrophy density Q(k,t) = nk(((k, t)((— k, £)), where 
(...) denotes averaging over thin circular shells: 

d t + 2uk 2 )n(k,t) =T n (k,t). (2) 

Here, the enstrophy transfer function Tq(A;, £) is given by 

T (M) = irk f P4^( C ( p ,t) C ( q ,t) C (-k,i))^ + c.c, (3) 

where c.c. stands for the complex conjugate term. Assuming that the system is in statistical 
steady state and extending integration in (H) only over all such triangles (k, p, q) that 



\k — p\ < q < k + p and p and/or q are greater than k c , one defines the two-parametric 
transfer T n (k\k c ) and TPEV 0: 

In a wide class of quasi-normal approximations the two-parametric transfer T^(k\k c ) 
in two dimensions is given by 



Ta(k\k c )= / / Q-k,p, q {p 2 - q 2 ) sina 



' ' ~ q2 -n( q )n(k) + fc2 ~f 2 n(pMfc) 



dpdq, (5) 



fc 2 g 2 ' /c 2 p 2 

where ®_ k ,p,q is the triad relaxation time. Here, the angle a is formed by the vectors p and 
q, and / J A denotes integration over the area defined above 

The main difference between various spectral closure models is in specification of Q-k, p , q - 
In [|I|], Tn(k\k c ) was evaluated using TFM. It was found that TPEV is a sign-changing 
function of the form u(k\k c ) = \u(0\k c )\N(k/k c ), with u{0\k c ) < 0, JV(0) = -1, and iV(l) w 
2.1. The derivation of 0-fc, P , g using the RG theory was given in |3[] and adapted for 2D 
isotropic and anisotropic turbulence in [|J and ||, respectively. In the present work, we 
compare u(k\k c ) for the inverse energy cascade regime calculated from DNS data with those 
predicted by TFM and the RG theory. 

We solve Eq. (0) numerically in a periodic box of the size 2n x 2tt using 512 2 Fourier 
modes. The numerical scheme involves a Fourier- Galerkin pseudo-spectral spatial approx- 
imation with implicit Adams- type second order stiffly stable time-stepping scheme 0. In 
order to increase the effective inertial range, mode-selective hyperviscosity of the form 
u(k) = fiik) + fs(k) = Aik~ 10 + Ask 14 has been introduced in the vorticity equation ([I]) 
instead of the molecular viscosity. The constant coefficients Ai and As have been selected 
empirically so as to minimize distortion of the energy inertial subrange. 

To simulate the inverse energy cascade, high wave number forcing 

f(k,t) = A f + 5 ktkf + 5 ktkf+1 ) e**® 
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is introduced in the vorticity equation; here, Af is the forcing amplitude, and random 
variable <f>(t) is a white-noise function of time t uniformly distributed in the interval [0, 2ir]. 
The results are not sensitive to initial vorticity field distributions, including the extreme 
case of the zero field. A series of numerical experiments with various kf and other flow 
parameters have been performed. Here we shall report only the results with kf = 98 since 
they gave the broadest inertial range. Other parameters used in these simulations were 
Af = 0.03244, A L = 0.5, and A s = 0.95 x 10" 34 . The value of A s chosen is somewhat high, 
to ensure efficient enstrophy dissipation and to eliminate the need for the dealiasing. 

In Figs la,b we plot the total energy E to t(t) = J +oc k~ 2 Q(k, t)dk and enstrophy VL to t(t) = 
Jo + °° Q(k, t)dk as functions of time. In Fig. la, one can see that the energy grows with time 
and eventually tends to reach a steady state. However, the drift towards the energy steady 
state is significantly slower than towards that of the enstrophy. Defining the rms velocity as 
Kms = l u (k)| 2 and the characteristic turnover time of the largest eddies as r tu = 2Tr/V Trils , 
we infer from Fig. lb that a steady state for the total enstrophy was achieved after about 
1.2 r tu , while about 5r tu were required to attain a steady state for the total energy. Note 
however that all the modes with k > 5 have reached the steady state after t « 2r iu , and 
only the largest modes were still developing. The results presented below pertain to the 
integration time t < 10 4 , before the energy saturates at low wave numbers. 

In Fig. 2 we plot the time-averaged energy spectrum obtained after about 5 r tu . The 
inertial range E oc k~ x extends over more than a decade in wave number space. Mean square 
line-fitting over the interval k e (12, 50) gives the scaling exponent close to the Kolmogorov 
value of |. Note that good agreement with the Kolmogorov scaling in the energy sub-range 
has been reported recently in || for 256 2 simulations and in || for very high resolution sim- 
ulations with 2048 2 Fourier modes. In Fig. 2 we also plot a compensated energy spectrum, 
k 5 / 3 e~ 2 / 3 E(k), where e is the energy transfer rate. The value of the Kolmogorov constant 
calculated from this data is about Ck = 6.2, in reasonable agreement with 5.8, calculated 
from DNS in |TIJ using the 256 2 resolution and 6.69, obtained analytically in |TT] on the 



basis of TFM. 



In Figs 3a,b we plot the fc- dependent energy and enstrophy flux functions, defined as 
n_B(fc) = Jq Tn(n)n~ 2 dn and Hn(k) = f Tn(n)dn, respectively As expected, an inverse 
energy cascade with constant energy transfer rate e develops for k < kf = 98, see Fig. 3a. 
For k > kf, TlE(k) and Hn(k) both quickly fall to zero, due to the strong dissipation at wave 
numbers adjacent to kf. In Fig. 3b, we see that the flux of enstrophy in the energy sub-range, 
k < kf, is zero. Strong enstrophy flux is observed for k > kf, until the enstrophy dissipation 
takes over and suppresses the flux of enstrophy into even smaller scales. The resolution 
employed in this study was insufficient to detect a well defined enstrophy transfer range. 
The results plotted in Figs 3a,b also indicate that the numerical scheme used conserves both 
total energy and enstrophy, since 11^(0) = 110(00) = and n B (0) = II^oo) = 0. 

We have calculated k and /c c -dependent enstrophy transfer function T^{k\k c ) employed in 
(d) by computing the third-order vorticity cumulants in (|3]) extending the integration only 
over those p and q that either p > k c or q > k c . We have set k c = 50, well inside the energy 
inertial subrange. 

The DNS-inferred normalized TPEV [ viz., the function N(k/k c ) = v(k\k c )/\u(0\k c )\] is 
plotted in Fig. 4, along with the TFM-based [TJ and RG-based || analytical predictions. 
The agreement between the DNS-based results and the TFM and RG theories is very good 
over the entire energy transfer range, up to the wave numbers close to k c , where the DNS 
data saturates, while TFM and RG curves exhibit sharp cusp. The physics leading to this 
cusp is as follows. As closer k approaches k c , as more elongated triads with either p or q <C k c 
become involved in the energy exchange between the mode k and the subgrid scale modes. 
The contribution from these triads results in the cusp behavior of the theoretical TPEV. 
However, in finite box DNS with large-scale energy removal, the energy of small wave number 
modes is reduced (see Fig. 2), which implies that instead of the sharp growth, the TPEV 
should saturate at k —>■ k c . To illustrate and quantify this explanation, we recalculated the 
RG-based TPEV with the enstrophy spectrum in (|5|) corrected at k < 5 according to Fig. 
2. In Fig. 5, we compare the DNS- and RG-based TPEV in their actual values, whereas the 
RG calculations were based upon the value of e found from DNS. The agreement between 
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the two is very good. We have also calculated TPEV for k c = 35,45 and 55 and found 
that the DNS-inferred TPEV scales with k~^ 3 , in full agreement with the Kolmogorov and 
Richardson laws. At all values of k c tested an equally good agreement between DNS data 
and RG predictions was observed. 

The good agreement demonstrated in Figs. 4 and 5 provides an indirect validation of 
TFM and RG results for isotropic 2D turbulence. 
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FIGURES 

FIG. 1. Evolution of the total energy E(k) (dotted line) and enstrophy Q(k) (solid line) towards 
the steady state. Dashed line denotes the total energy with the first 6 modes excluded. 

FIG. 2. Energy spectrum E(k) (solid line) and compensated energy spectrum E(k)k 5 / 3 e~ 2 / 3 
(dotted line). 

FIG. 3. The energy flux n#(fc) (solid line) and the enstrophy flux ILn(k) (dotted line). 

FIG. 4. Normalized two-parametric eddy viscosity from DNS (dots), from TFM (dashed line), 
and from RG (solid line). 

FIG. 5. Actual two-parametric eddy viscosity from DNS (dots) and from RG (solid line). In 
RG calculations, the energy spectrum for k < 5 was corrected in accordance with the DNS results, 
Fig. 2. 
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